Universal discriminator for completely unknown optical qubits 
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We propose an experimental setup that is capable of unambiguously discriminating any pair of 
linearly independent single photon polarization qubits, about which we don't have any knowledge 
except that an extra pair of these unknown states are provided as the reference. This setup, which 
is constructed with optical CNOT gates, weak cross Kerr non-linearities, Bell state analysers and 
other linear optical elements, transforms the unknown triple photon input states to the correspond- 
ing single photon states to be deterministically processed by linear optics circuit. The optimal 
discrimination of the unknown states is achieved by this setup. 
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I. INTRODUCTION 

When we discriminate a pair of completely unknown 
quantum states for i = 1,2, the only possible way 
is to have copies of these unknown states as the refer- 
ence, and then we can design the proper measurement 
that is capable of distinguishing between the inputs, e.g., 
|V'i)IV'2)|V'i)> in which the data states li/^i) are appended 
to the reference copies to form the states with the per- 
mutation symmetry. If we require that the incoming 
unknown data states should be unambiguously discrim- 
inated, our device will output three possible results: 1 
corresponding to the input being iV'i), 2 to that being 
\ip2) , and for inconclusiveness. The optimal unam- 
biguous measurement of these symmetric input states has 
been studied in both the Bayesian [1] and the minimax 
approaches 

Another feature we need to have for such a device is 
its universality: it will perform optimally for any pair of 
unknown states. Since the input states are unknown, the 
only quantity we can use to indicate how well the device 
works is the average success probability over the Hilbert 
spaces of all randomly distributed lipi) and \tp2), which 
are discriminated by this single device. For the discrim- 
ination of two unknown qubits with one copy of them as 
the reference, the optimal average succes probaility can 
be as large as 1/6, if these two unknown qubits appear 
with an equal a priori probability [l[ . 

There are ways to improve this optimal average success 
probability. We can prepare multiple copies of reference 
and apply different types of measurement, i.e., projection 
value measure (PVM) and positive operator value mea- 
sure (POVM) working in different range of a priori prob- 
ability of the inputs, to significantly increase this value, 
if we follow the Bayesian approach to deal with the prob- 
lem [Sj] . For the discriminators to unambiguously identify 
the averaged unknown input states, the IDP (Ivanovic- 
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Dieks- Peres) upper bound [J, la, |6[ of success probability 
in discriminating a pair of known states can be achieved 
if the copy number of the reference states goes to infinity 
[j]. To increase the optimal average success probabil- 
ity further, we can use both multiple reference and data 
copies, so that the device will realize the deterministic 
discrimiantion of any pair of unknown states in the limit 
of infinite reference and data copy numbers together 

This type of unknown state discriminators can be of 
potential value in quantum communication and quantum 
computing, so it is interesting for us to study their fea- 
sible physical implementations. In [9| we propose a pure 
optical way based on interferometry [lO| to deterministi- 
cally realize the optimal measurements in IH under the 
assumption that all the symmetric input states can be 
prepared as single photon states in multi-rail representa- 
tion. The quantum circuit to realize the unknown qubits 
discrimination without achieving the optimal measure- 
ment and the possible implementation of its ingredient 
logic gates in ion traps are studied in [ll[ . The design of 
the unitary operators to implement the discrimination of 
a pair of unknown qubits is also discussed in [T^ . 

In this paper, we design an experimental setup to un- 
ambiguously and optimally discriminate the symmetric 
input states jfAi) IV'2) iV'i) (* = 1:2) prepared with any 
pair of linearly independent unknown single photon po- 
larization qubits. This is a device to process the unknown 
triple photon input signals, and all the technologies in- 
volved have been experimentally realized thus far. After 
the transformation from the triple photon inputs to sin- 
gle photon signals, only the ancilla of vacum state will 
be used in the realization of the optimal POVM the de- 
vice performs. The rest of the paper is organized as fol- 
lows: in Sec. II we present the physical basis to realize 
the optimal measurements designed to discriminate the 
general unknown qubits, which we intend to realize for 
the polarization photon systems; in Sec. Ill we look at 
the state preparation part of our device, where the triple 
photon input states are transformed to the single photon 
states through a combined teleportation protocol, so that 
they can be deterministically processed by multi-port in- 
terferometer; the simplest linear optics circuit to finally 
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determine the incoming unknown date qubit is given in 
Sec. IV; and finally we give some conclusive discussions 
in the last section. 

II. PHYSICAL BASIS TO DISCRIMINATE 
UNKNOWN OPTICAL QUBITS 

We first give a brief review of the optimal measure- 
ments for discriminating any pair of unknown states we 
aim to realize for single photon polarization qubits. The 
input signals processed by our device are given as follows: 

\'^l,^n) = |V'l)l|V'2)2|V'l)3 

|^2,») = |^l)l|V'2)2|V'2)3, (2.1) 

where {ipi) = oii\H) + l3i\V), for i — 1,2, and H, V are 
their horizontal and perpendicular polarization modes 
[l3l |. Since j^/'i) and |'!/'2) are linearly independent, there 
exist different permutation symmetries in these two in- 
puts, and we can therefore design the POVMs, with the 
elements Hi and U2 satisfying Hil^-I") = nal^T) = 0, 
to unambiguously discriminate them. In the Bayesian 
approach [l|, the optimal POVM elements for the unam- 
biguous discrimination of the above inputs are 

= ^(^2-y|)/2®|vI'-)(vI/-|, (2.2) 



where l^f^) = l/V2{\H),\V)j-\V),\H),) and L the 
identity operators, if the a priori probability rji jl4l| of 
the data \tpi) is in the range of 1/5 < r/i < 4/5. For 
4/5 < f^i < I7 the optimal measurement reduces to the 
PVM as 

Hi = mi) (nil 

II2 = 0; (2.3) 

and for < 771 < 4/5, the index 1 and 2 of the above 
PVM operators should be interchanged to get the op- 
timal measurement. In the minimax approach without 
the available a priori information of the input data 0] , 
on the other hand, the optimal POVM for discriminating 
the inputs in Eq. (2.1) is given as 

n2 = ^hf^MDimi (2-4) 

In both approaches, there is inconclusive measurement 
result corresponding to the operator Hq = / — Hi — 112. 

The physical fundamentals to realize the optimal mea- 
surements on the unknown input signals are the telepor- 
tation of unknown quantum states [15] and the determin- 
istic implementation of a POVM on single photon signals 
[T^ . A teleportation protocol is applied to transform the 
triple photon input signals in Eq. (2.1) to the correpond- 
ing single photon signals running on different tracks: 



|*l,m) = |^l>l|^2>2|^l>3 

ci\HHH) + C2\VHH) + cslHHV) + ■■■ + cglVVV) 
^ ci\0a,0b,H)+c2\Ia,0b,H)+c3\0a,0b,V) + --- + cs\1a,Ib,V), (2.5) 



where Ci are unknown coefficients, and Oa, 1a, Ob and 
1 B the symbols related to the which-path degree of free- 
dom of the photon. After some of the input signal's po- 
larization degree of freedom is mapped this way to the 
which-path degree of freedom, we will obtain a single po- 
larization photon running on 4 different tracks (the total 
dimensionality of the quantum system is 8) so that any 
POVM performed on it can be deterministically imple- 
mented by linear optics. 

To teleport three input qubits in a combined way, we 
need to have three entangled pairs, e.g., 

\^t) = -^i\H)\H) + \V)\V)), 

= -^i\H)\OA) + \V)\lA)), 

1$+) ^ J-(\H)\Os) + \V)\Ib)). (2.6) 



We here define the four Bell states as follows: 
|<i>±) ^ i=(|0)|0)±|l)|l)), 

|vl.±) ^ i=(|0)|l)±|l)|0)), (2.7) 

where can be H, Oa, Ob, etc, and 1 can be V, 1a, 
1b, etc. Two of the states in Eq. (2.6) involve both 
polarization and which-path degrees of freedom. How to 
effectively produce them is what we will discuss in the 
next section. 



III. CONVERSION TO SINGLE PHOTON 
SIGNAL 

In this section, we give a rather detailed description of 
a circuit to map the triple photon input signals to the 
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corresponding single photon signals as in Eq. (2.5). 




To begin with, we look at the transformation of a dou- 
ble photon polarization state \^in) = |V'i)i|V'2)2, where 
lipi) — ai\H) + (3i\V), to the corresponding single pho- 
ton state by a combined teleportation procedure. We 
use two extra photon sources 1/^/2{\H)a + \V)a) and 
1/V2{\H)b\H)c + \V)b\V)c), and let the first photon 
on track A control the part B of the second entangled 
photon pair through a CNOT gate taking the following 
action: 



FIG. 1: (Color online) The quantum non-demolition detec- 
tion to seperate even parity {\H)b\H)c,\V)b\V)c} from the 
odd parity {1-^)51 V)c|V^)_b| H) c } states. Through one of the 
weak non-linear cross-Kerr, the laser probe initially in the co- 
herent states |a) will evolve to |ae'^) or |ae~'^), if the single 
photon is present too. By a homodyne detection |ae*'^) re- 
sulting from the odd parity can be discriminated from |a) 
from the even. We classically feedforward the detection re- 
sult to the switch, S, to direct the odd parity component to 
the track Bi and the even parity component to the track B2 . 
On Bi the odd parity is restored to the even by the bit flip 
and the necessary phase shifters. We use a fictitious control- 
ling action qubit, a\OA} + I3\1a), acting on the track B photon 
component to describe the combined action of the CNOT gate 
and this QND. 



J 



UcNOT = {\H)a{H\ <E>Ib + \V)a{V\ <E) a^,B) ^ Ic- (3.1) 



The above equation means that there will be no action 
on track B if the photon component on track A is in the 
state of \H) and there will be a bit flip ax.B of {H,V} 
on track B if that photon component is in \ V) . After the 
joint unitary map Ucnot, the total state on the track A, 
B and C will be as follows (here we neglect the common 
factors for brevity): 



Ucnot{\H)a + \V)a){\H)b\H)c + \V)b\V)c) 

\H)a\H)b\H)c + \H)a\V)b\V)c + \V)a\V)b\H)c + \V) a\H) b\V) c ■ 



(3.2) 



r 



This type of non-destructive optical CNOT gate has been 
experimentally demonstrated [13] and, theoretically, a 
CNOT gate for two independent single photon polariza- 
tion qubits can be realized near deterministically [18]. 

The output state of the CNOT gate in Eq. (3.2) 
is divided into two sets for part B and C: the 
even parity {\H) b\H) c ^\V) b\V) c} and the odd parity 
{\H)b\V)c,\V)b\H)c}, and they can be discriminated 
near deterministically by a polarization parity quantum 
non-demolition detection (QND) [l^. The setup to real- 
ize this QND includes a laser probe initially in a coherent 
state \oi)p and the cross-Kerr nonlinearities, which have a 
Hamiltonian H = hxalaga^ap, where the signal (probe) 
mode has the creation and destruction operators given 
by a|, as{ajj, ap) and x the strength of nonlincarity which 
can be weak. Through the interaction in the weak non- 
linear cross-Kerr in Fig. 1, the photon states on track B, 



C and the probe evolve together to ]'0t) ~ i\H)B\H)c + 
\V)B\V)c)\a)p + \H)B\V)c\ae'')p + \V)B\H)c\ae-'')p, 
where 9 = xt with t being the interaction time. \a) and 
]Q;e^*^) can be distinguished between each other by a 
homodyne-heterodyne measurement on the probe. 

Then we classically feedforward two actions to the out- 
put on track B: if the the detection result is ]Q!e^*^), 
we let the photon component on track B go to path Bi 
and denote the action 0^ (the index A means the con- 
trol from track A photon component); and if the result 
is ]q!), the photon component on track B is redirected 
to track B2 with an action named I a- A bit flip (plus 
the phase shift redressing due to the QND) on Bi track 
is performed to restore the odd parity of the photon 
components on track B and C to the even parity, i.e., 
\HV), \VH) \HH), \VV). As the result, we obtain a 
state, 
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1^^"*) ^ \H)A\H)B,flAH)c + \H)a\V)b,.oAV)c + \V)a\H)b.,iAH)c + J^)c 

= mA\QA) + |l^)A|U))(|i/)B|H)c + \V)b\V)c). 

I 



(3.3) 



which can be factorized into the product of two entan- 
gled states. \Qa) and are defined as two controlling 
switch actions on the track B photon component from 
the photon component on track A. These two fictitious 
quantum states describe the effective action of track B 
photon component to go to different paths Bi and B2 
under the control: 

\Qa)\H)b ^ \H)b,,oa, \W\H)b ^ \H)b,,1a, 



\0a)\V)b^\V)b,,Oa, \Ia)\V)b^\V) 



(3.4) 



It is feasible, therefore, to use these two entangle states 
to teleport l^"™) = l'02)2 to a single photon polar- 
ization state: 



\i^l)l{\H)A\OA) + \V)a\1a))H2)2{\H)b\H)c + \V)b\V)c) 

(|$+)l,A|V'l)Act + I*+)i,a(CT:.,AcJ^1>AcJ + \'^')l,AH(^y.AcMl)Act) + \'^')l,Ai<^z,Act\i^l)Act)) 
X (|$ + )2,c|^2)b + |* + )2,c(CT:.,i3|V2)s) + | 2,C (^io-y ,B | V'2> S ) + | ) 2,C (^^.S | V'2> S )) , (3.5) 

I 



where the controlling action qubit is defined as \tpi)Act = 
ailOyi) + PiI^a), and ax, Cy and Oz the Pauli matrices. 
With a final restoring unitary transformation Urest on 
both the tracks Si, B2 and the polarization modes H , V 
according to the information fed forwarded after the Bell 



state analysis on track 1, A and 2, C, respectively, we 
will obtain a single photon polarization state running on 
track Bi and B2 (the total dimensionality of the outputs 
is 4 as the inputs): 



Infant) = \^l)AcAlb2) ^ {ai\OA) + Pl\lA)){a2\H)B + f32\V)B) 

= aia2\H)Bi,0A + aiP2\V)Bi,aA + l3ia2\H) B2.IA + l3il32\V) b2,i. 



(3.6) 



r 



r 



The complete (and close to deterministic) Bell state 
analysis setups can be found in [l^, [2l| , and also in 
[22,] a practical method for teleporting optical qubits is 
given. These experimental schemes can be modified to 
realize the operations in our setup. The whole scheme to 
realize the transformation from double to single photon 
states is given in Fig. 2. 

So far we have put together all necessary ingredients 
of a circuit to transform the triple photon inputs in Eq. 



(2.1) to the corresponding single photon states. To tele- 
port the triple photon inputs to the corresponding single 
photon states, we just need to continue the above proce- 
dure before the Bell state analysis by adding one more 
single photon source, 1/\/2{\H)d + \V)d), and two more 
CNOT gates and parity QNDs. We use a 50-50 beam 
splitter to split the second single photon into two differ- 
ent paths: 



\H) 



\V)i 



\H)d, + \V)d, + \H)d2 + \V)d2, (3.7) 
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FIG. 2: (Color online) The scheme to transform a two-photon state, I^Pin) ~ \il^i)\ip2), to a single photon polarization state 
running on two different tracks, I'ifout)- We use a BBO to generate an entangled pair l/^/2{\H) b\H) c + \V)b\V)c) on track 
B and C, and put a single photon 1/V2{\H)a + \V)a) on track A. The CNOT gate can be a modified version of the cited, and 
the QND unit is shown in Fig. 1. The results of two Bell state analysis are sent through the classical channels to where the 
restoration unitary map Urest is performed. Urest is chosen from a finite set of operations. 
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FIG. 3: (Color online) The scheme to teleport any symmetric unknown three-photon input state |i/'i)|i/'2)|i/'i), for i = 1,2, 
to the corresponding single photon state. One more single photon source is used on track D, and it is split by a 50-50 beam 
splitter BSi (and the necessary phase shifters) into two paths to control the respective photon components. Another identical 
beam splitter BS2 merges the photon components on track Di and D2 to a single path again. 



and the photon components running on Di (-D2) track 
controls the photon component on Bi (-B2) track through 
an optical CNOT gate. Then, performing the polariza- 
tion parity QNDs on track Bi, C and B2, C, respectively, 
and the necessary bit flips and phase redressing, we will 



seperate the photon components on Bi and B2 to those 
on 4 tracks -Bi.i, -61,2, B2.1 and -62,2, and thus obtain the 
output state which can be factorized into the product of 
three entangled states (similar to the state in Eq. (3.3)). 
We are now able to teleport 3 unknown qubits together 
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to a single photon polarization state running on 4 differ- discriminate the data \ipi) and \tp2) 
ent tracks. The complete scheme of the teleportation of 
the triple-qubit to a single photon state is shown in Fig. 
3. 



IV. OPTIMAL POVM IMPLEMENTATION 
CIRCUIT 



The circuit in Fig. 3 generates the following single 
photon states as the inputs of the POVM modules to 



= {ailOo) + Pi\lD)){a2\0 a) + P2\l a)) {a^\H)B + f3i\V)B) 

= aia2ai\0D,0A,HB) + aia2/3z|0_D, 0^, Vb) + ail32ai\0D,lA, Hb) + Q;i/32/3i|0D, 1a, Vb) 

+ f3la2a^\lD,0A,HB) + f3ia2Pi\lD,0A,VB) + Pi(32a,\lD,lA,HB) + PlP2P^\lD,lA,VB) 

= aia2ai\H)Bi^i + aia2Pi\V) Bi,i + aiP2ai\H) b2a + aiP2Pi\V) B2.1 

+ /3ia2a,\H)B,,2 + Pla2/3^\V)B,,2 + ^3lP2a^\H)B2,2 + /5i/52A|V")i3.,2, (4.1) 



Mirror 




(N-iy 



FIG. 4: (Color online) The unitary transformation mod- 
ule constructed with beam splitters (dark square) and phase 
shifters (white square). Any unitary operator represented by 
an N X N matrix can be implemented by an TV x A*' ports 
module of this kind with the maximum N{N — l)/2 of beam 
splitters p^ . 



where the actions defined in Eq. (3.4) are used in the 
notation. With 4 polarization beam splitters (PBS) we 
can thus obtain the multi-rail single photon signal states 
running on 8 tracks. Any POVM with n elements n^, 
for i = 1, • • • , n, on such single photon states can be 
deterministically implemented with at most 3n — 2 linear 
optics modules shown in Fig. 4 [l6l |. 

The POVM for unambiguously discriminating these 
states satisfies more condition, Hij^'j j„) = 112 1^*! in) = 
0, so we can use a linear optics scheme, which requires 
fewer and much simpler modules, to realize the unam- 
biguous discrimination of and 
scheme works with the following unitary map Q 



Such a 



Un 



odule 



(/-no) 



(4.2) 



which is realizable by one 16 x 16-port module shown in 
Fig. 4. This module transforms the signals embeded in a 
larger space of the double of the signal space dimension, 
(l^i j„), 0)"^, to the components in two orthogonal output 
subspaces: 



Ur, 



odule 



I*' 



i.inl 





(4.3) 



The component (/ — ^o)'^\^'iin) corresponds to the 
success of the measurement with a probability P = 

{K^n\I - no|*^,.„) = {%^n\^^\%^n), and the COm- 
1 

ponent ngl^'Jin) to the inconclusive result, which oc- 
curs with a probability Q — (^i,i„|no|^- j„). It is easy 
to verify that |5''ij„) and 1^*2 i„) are mapped to distin- 
gushable orthogonal states in the first subspace, since 

.> = o. 



2,inl 



The linear optics circuit to perform Umoduie can be 
simplified further. We first decompose the inputs in 
Eq. (4.1) into the components in the orthogonal sub- 
spaces Hi = {\0d,0a,Vb),\0d,1a,Hb)A^d,0a,Hb)}, 

n2 = {\lD,OA,VB),\OD,lA,VB},\'i^D,U,HB}} = 

{\Qd,Oa,Hb)} and H4 = {\IdAa,Vb)}, respectively. 
We introduce a unitary map Vi, 
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|*i> = ^{\0d,0a,Vb)~\0d,Ia,Hb)) 

\'^2) = ^{\0d,0a,Vb) + \0d,Ia,Hb)) ~ ^\1d,0a,Hb) 

1*3) = ^J^{\0d,0a,Vb) + \Od,Ia:Hb) + \Id,Oa,Hb)), 
in Hi, and another unitary map V2, 

= ^JI{\1d,0a,Vb)-\1d,Ia,Hb)) 

1*2) = ^JI{\Id,Oa,Vb) + \1d,Ia,Hb)) - ^\1d,0a,Vb) 

\%) = \II{\Id,0a,Vb) + \1d,Ia,Hb) + \0d,Ia,Vb)) 



(4.4) 



(4.5) 



r 



in 7^2, respectively. These two unitary maps on the 
modes of the relevant tracks correspond to two simple 
and identical 3 x 3-port modules processing the rele- 
vant components of the input signals. In the subspaces 
n[ = {|$i), 1*2)} and n'^ = 1$^)}, Ho of the op- 

timal POVM is given as a 2 x 2 matrix [1]: 



n, 



(■0 



(1 



(2- 




for i = 1,2, and the total Hq is the direct sum 



Hn = n, 



(1) 



(2) 



(4.6) 



(4.7) 



Therefore, the signal components outside H'l and 7^2 only 
contribute to the inconclusive results, and we only need 
to process the signal components in these two subspaces. 

For the a priori probabilities r/i = r/2 = 1/2, i.e., the 
most difficult situation to discriminate the input data, 



as well as the optimal POVM obtained in the minimax 

approach the operator IIq^'' and IIq^'' are therefore 
reduced to 



n, 



(1) 



n 



(2) _ 



2 

h 

6 



6 

2 
3 



(4., 



The unitary transformations in Eq. (4.2) constructed by 
these two non-unitary operators are, respectively, imple- 
mented by two identical circuits consisting of only 3 beam 
splitters, one of which performs the unitary map. 




(4.9) 



which diagonalizes IIq^'' and TI^q \ and the other two of 
which perform the unitary map in the extended space, 




(4.10) 



Through all these unitary transformations the different spectively, have been mapped to the orthogonal states, 
input signal components from l^'i j„) and |5'2 m)j re- 
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FIG. 5: (Color online) The layout of the optimal POVM im- 
plementation. The input |^i„) is the output \'^out) in Fig. 
3. The polarization beam splitters (PBS) seperate the input 
into the 8 dimensional single photon state in multi-rail repre- 
sentation. The Vi are 3 x 3-port modules, and Ui i x 4-port 
modules with two input ports black. Only 4 photon detec- 
tors are required to place on the terminals where the useful 
information is output. At all other terminals (2 of the PBS, 
2 of Vi modules and 4 of Ui modules), there are only useless 
signals contributing to the inconclusiveness. We can set the 
output in such a pattern: if the detector Ci or C3 clicks, we 
conclude that the unknown input is and if C2 or C4 

clicks, the input date must be IV'2). From each counter, we 
detect the photon with a probability 1/12, if the data 
and \tp2) occur with the same a priori probability. 



Adding one more 50-50 beam splitter to seperate the sig- 
nals from the different inputs, we realize the target to 
discriminate the data input: the input data is determined 
to be l-^i) if the output photon is recorded by one single 
photon detector and \ip2) if the outout photonis recorded 
by the other. In Fig. 5 we use Ui to represent a whole 
package of these unitary maps including Ui^i and Ui^2- 

By a straightforward calculation considering Eq. (4.4)- 
(4.5), (4.9) and (4.10), we see that in one of the photon 
detectors in Fig. 5 there is an average probability of 1/12 
to detect the photon. With the detection of the photon, 
we can draw the conclusion about which the input data 
is. Adding together the contribution of another photon 
detector placed at one terminal of the other Ui module, 
we obtain the optimal average success probability 1/6 in 



the discrimnation of a pair of unknown qubits [l[. The 
layout of the optimal POVM implementation is given in 
Fig. 5. 



V. SUMMARY AND DISCUSSION 

We have constructed a setup to unambiguously and 
optimally distinguish between any pair of linearly in- 
dependent single photon polarization qubits, which are 
completely unknown to us. This proposed experimen- 
tal setup applies the state-of-the-art technologies to pro- 
cess quantum information, e.g., the non-destructive op- 
tical CNOT gates, the complete and near deterministic 
Bell states analysis including the teleportation of opti- 
cal qubit, etc.. which have been experimentally realized 
0, E 0, ml HI- Neglecting the possible error in the 
operation, the optimal discrimination of any pair of un- 
known optical qubits can be deterministically performed 
with this device. The future advance in these techonolo- 
gies will make such a device more reliable and less costly. 

In principle, we can also generalize the setup to a de- 
vice that works with multiple reference or data copies, 
which will increase the success probability in discrimi- 
nating the unknown data 0], or a device that discrim- 
inates more than two unknowm states This can 
be done by teleporting more controlling switch action 
states^ but the restriction is the cost of the above men- 
tioned experimental facilities. Another generalization is 
a setup that unambiguously and optimally discriminates 
any pair of unknown d-level quantum systems, i.e., qu- 
dits. Due to the lack of knowledge about the quantum 
circuits of the qudits, it seems that so far there has been 
no available experimental facility to directly discriminate 
two unknown qudits unless we first transform them, to- 
gether with the reference copies, to qubits with a certain 
possibility. These interesting problems, as well as the 
experimental implementation of the setup we propose in 
this paper, are worth further study. 
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